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ABSTRACT 

At present, 30-40 per cent of the baryons in the local Universe is still undetected. Accord- 
ing to theoretical predictions, this gas should reside in filaments filling the large-scale struc- 
ture (LSS) in the form of a Warm-Hot Intergalactic Medium (WHIM), at temperatures of 
10^ — lO*^ K, thus emitting in the soft X-ray energies via free-free interaction and line emis- 
sion from heavy elements. In this work we characterize the properties of the X-ray emission 
of the WHIM, and the LSS in general, focusing on the influence of different physical mecha- 
nisms, namely galactic winds (GWs), black-hole feedback and star-formation, and providing 
estimates of possible observational constraints. To this purpose we use a set of cosmologi- 
cal hydrodynamical simulations that include a self-consistent treatment of star-formation and 
chemical enrichment of the intergalactic medium, that allows us to follow the evolution of 
different metal species. We construct a set of simulated light-cones to make predictions of the 
emission in the 0.3-10 keV energy range. We obtain that GWs increase by a factor of 2 the 
emission of both galaxy clusters and WHIM. The amount of oxygen at average temperature 
and, consequently, the amount of expected bright Ovil and Ovill lines is increased by a factor 
of 3 due to GWs and by 20 per cent when assuming a top-heavy IMF. We compare our results 
with current observational constraints and find that the emission from faint groups and WHIM 
should account from half to all of the unresolved X-ray background in the 1-2 keV band. 

Key words: Cosmology: theory, large-scale structure of Universe - X-rays: diffuse back- 
ground, galaxies: clusters - methods: hydrodynamical simulations 



1 INTRODUCTION 

The amount of baryonic matter present in the Universe is nowa- 
days measured with high accuracy in the framework of the stan- 
dard cosmological A-cold dark matter (AC DM) scenario. Accord- 
ing to the most recent CMB observations jKomatsu et alJ[2oTlh 
the percentage of baryonic density with respect to the critical one 
is 100 r^i, = 4.49 ± 0.28, a figure which is in excellent agree- 
ment with the estimates comi ng from the primord ial nucleosyn- 
thesis of the heavy elements jKirkman et al.|[2003h and with the 
baryo n budget at high redshift estimated from Lya absorption sys- 
tems ^Weinberg et aI.|[l997l ; lRauchlll998l) . thus representing one of 



the most important result of modem cosmology. Even if they con- 
stitute a minor component in the matter-energy budget, baryons are 
indeed crucial because, since they interact with the electromagnetic 
field, they are responsible of all the radiative phenomena, thus pro- 
viding us the only probe to unveil also the dark components. 

On the other side baryons associated to Lya absorbers at 
the pr esent epoch account for only ^ 30 per cent of the cosmic 
mean jPenton et al.ll2004 ISregmanllzOOTb . Indeed, when the cos- 
mic structures form, the baryon census becomes much more chal- 
lenging and it requires observations at different wavelengths. De- 
spite the great improvements in our observational capabilities of 
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the last decade, when accounting for all the matter observed from 
the optical band (stars) to the X-rays (clusters) at z < 2 the amount 
of observed gas in the Universe is much lower than expected. More 
precisely, about 30-40 per cent of the gas is still out of reach from 
current observational instruments (see lNicastro et al]|2005al for a 
review): this issue is currently referred to as the "missing baryons" 
problem. 

The first attempts to address this matter from the the- 
oretical point view come from hydr odynamical simulations 
dCen & Ostrikerlll999l : iDave et al.|[200lh that showed that when 
galaxy clusters form, a significant amount of the whole intergalac- 
tic medium (IGM) remains out of the virialised halos creating a 
network of filaments, the so-called cosmic web ( lBondetal.ll99d) . 
that provide a continuous flow of accreting material to the clusters 
themselves, located at the knots of the structure. The gas in these 
filaments remains at moderate overdensities (5 ~ 10 — 10'^) and 
is shock-heated up to temperatures of the order of 10^ — lO'^ K, 
and is usually addressed as the Warm-Hot Intergalactic Medium 
(WHIM). Given its thermodynamical conditions the hydrogen and 
helium are fully ionized and emit via free-free interaction with the 
electrons. This process together with the line-emission from partly 
ionized metals is expected to contribute to the X-ray background 
(XRB) in the form of a soft diffuse component. 

Analysing the current observational constraints, the unre- 
sol ved emission in the Chandr a Deep Fields (CDFs) observed 
by iHickox & MarkevitcH ( l2007h in the 0.5-1 keV band is con- 
sistent with the expectations o f the diffuse gas obtain ed from nu- 
merical simulations (see, e.g., iRoncarelli et"al]|2006l) . Moreover, 
lGaleazzi"et aP (l2009l) detected the signature of a possible WHIM 
emission in the correlation function of the diffuse emission signal 
in the 0.4-0.6 keV band obtained from several XMM-Newton obser- 
vations. On the other side, there have also been claims of detection 
of WHIM emission associated to local overdensities of the LSS. 
IWemer et aP ( l2008h observed at 5a a continuum emission from a 
filament extending between the Abell 222 and Abell 225 clusters. 
Similarly, ROSAT observations allowed the discovery of a soft dif- 
fuse emission consistent wi th T 10^ K gas ass ociated with a 
galaxy group at z ^ 0.45 jZappacosta et al.ll2003) and with the 
Sculptor supercluster jZappacosta et al.ll2005l) . 

Probably the most important point in WHIM observation 
strategies is the fact that this medium is believed to contain met- 
als whose emission and absorption lines provide an effective way 
to disentangle its signal from background and foreground objects. 
In particular, the best targets are C, O, Ne and possibly also N, Mg 
and Fe, with bright features expected from K and L-shell transi- 
tions of H-like, He-like and Li-like oxygen ions, respectively the 
OVIII Is - 2p doublet (E =654 eV), the Ovil Is — 2p resonance 
line (574 eV) and the Ovi doublet (12 eV). This method already 
provided several claimed detec tions, all regarding absorption lines. 
iDanforth & ShuH ( I2005ll200^ detected 83 Ovi absorption systems 
in the UV spectra of 43 AGNs. However, since these absorbers are 
associated to Lya counterparts it is no t clear how much of this gas 
can be a ccounted to the miss ing part jTripp et al]|2008h . At X-ray 
energies ' Nicastro et al.l l l2005ta' ) detected several OviI absorption 
features in the spectra of the Mrk 421 blazar, but the statistical sig- 
nificance of their measuremen t s has been quest ioned JKaastra et al.l 
l2006l : iRasmussen et alj 12007: Yao et al. 2011). Several other de- 
tections with both Cha ndra (Fang et al. 2002 ; Mathur et al. 20031 : 
Zappacost a et al.' '2 010l) and XMM-Newton (lFu|imoto et alj |2004| : 
Takei etalji2007. : ,Buote et all l2009l : iFang et all 1201 Oh have been 
claimed but their statistical significance is low and they mainly cor- 
respond to the densest part of the IGM (galaxy clusters or groups), 



indicating that these measurements are at the limit of current obser- 
vational capabilities. 

Up to now the detection of emission lines has proven to be 
even more difficult. However, detecting the WHIM in emission can 
be much more rewarding since it has the potential to map large 
structures and does not depend on the presence of background 
sources, thus potentially providing a global picture of WHIM 
thermo- and chemo-dynamics, which is by now far from being 
reached. 

In this scenario a big effort has been put into the modelisa- 
tion of the WHIM which becomes crucial in the perspective of fu- 
ture X-ray missions in particular for what concerns the treatment 
of radiative cooling, stellar and galactic feedback, reion ization and, 
most importantly, c hemical enrichment. lCen & Ostriker. (2006) and 
ICen & Fandtoeh included for the first time the enrichment from 
Type II supernovae (SN-II) of different metal species in their Eu- 
lerian cosmological simulations. Their work showed how the ef- 
fect of different galactic wind models can change both the ther- 
modynamics and_the abundance outside clusters. More recently 
iBertone et al. I l20l3) used a similar approach with a set of sev- 
eral large-scale simulations performed with the GADGET- 3 code, 
including also Type la supernovae (SN-Ia) and asymptotic giant 
branch (AGB) stars to study specifically the properties of metal line 
emission. 

In this w ork we use the set of c osmological hydrodynamical 
simulations of'Tomatore et al.' ('2010') to study and characterize the 
emission of the WHIM, with particular emphasis on the of impact 
different enrichment schemes. Starting from the outputs of these 
simulations, we adopt a light-cone reconstruction approach in or- 
der to provide a more direct comparison with current and future 
observational constraints. 

This paper is organized as follows. In the next s ection we will 
describe the simulation set of iTomatore et alj ( |2010|) which consti- 
tutes the basis of our study. In Section|3]we will explain our method 
to create the mock light-cones and spectra from the outputs of the 
hydrodynamical simulations. Section |4] is devoted to the analysis 
of the global properties of the IGM and WHIM emission, while 
Section |5] focuses on our predictions for the oxygen line statistics. 
We summarize our results and draw our conclusions in Section|6] 
Unless stated otherwise, when referring to solar metal abundances 
and relative element yields we will assume the results reported in 
lAsplund et al] ( l2009h . 



2 THE HYDRODYNAMICAL SIMULATIONS 

In this section we review the simulation set of iTornatore et al.l 
(2010). We refer the reader to their original work for further de- 
tails. The set consists of 8 cosmological hydrodynamical sim- 
ulations performed with the Tree PM- smoothed p article hydro- 
dynamics (SPH) code GADGET-2 dSpringell [2005h in a version 
that includes the implementation of the chemical enrichment of 
Tornatore et al. ( 2007). A summary of the main characteristics of 
the different runs is shown in Table [T] 

For all the runs the cosmology assumed is a flat standard 
ACD M model consistent with the WMAP-5 results dKomatsu et al.l 
|2009|) . namely with f7A=0.76, 0^=0.24 and r2b=0.0413 for the 
matter-energy density parameter of the cosmological constant, to- 
tal matter (dark and baryonic) and baryonic matter alone, respec- 
tively. The other cosmological parameters are a primordial power 
spectrum of the cold dark matter density perturbations with slope 
Us = 0.95 and normalization corresponding to (Tg=0.8 and a Hub- 
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Table 1. Characteristics of the 8 simulations object of our study. Column 1: run name. Column 2: comoving length of the box (h~^ Mpc). Column 3: number 
of particles. Column 4: number of available snapshots in the interval < 2 < 1.5. Column 5: stellar IMF assumed. Column 6; feedback scheme (see more 
details in the text). 



Run name 


^box 




Snapshots 


IMF 


Feedback model 


W (reference) 


37.5 


256^ 


25 


Kroupa 


Winds 


W75,512 


75.0 


512^ 


6 


Kroupa 


Winds 


W37_400 


37.5 


400^ 


16 


Kroupa 


Winds 


Way 


37.5 


2563 


25 


Arimoto-Yoshii 


Winds 


Ws 


37.5 


2563 


25 


Salpeter 


Winds 


NW 


37.5 


2563 


25 


Kroupa 


No winds 


cw 


37.5 


2563 


25 


Kroupa 


Coupled winds 


BH 


37.5 


2563 


25 


Kroupa 


BH-feedback, no winds 



ble constant Ho =100 h km s"^ Mpc"\ with h = 0.73. The grav- 
itational softening is set to e = 7.5h^^ kpc in comoving units at 
z > 2 while at lower redshifts it is set to e = 2.5h~^ physical kpc. 

We take the W run (see Table [T} as a reference model for 
comparison. This simulation follows the evolution of a box of 
L =37.5 comoving Mpc per side, with 2x2563 particles for 
dark matter and gas: this results in a mass resolution of mdm = 
1.9xl0®fe"^ M0 andrrigas = 3.6x lO^/i"^ Mq, respectively. The 
choice of these parameters is the result of a compromise between 
the need to sample the perturbations on large scales and to achieve 
the resolution necessary to describe the dense environments where 
star formation and feedback take place. This run produced 43 out- 
puts logarithmically equispaced in redshift from 2 = 4 to z = 0. 
Since we are interested mainly in the WHIM emission and we do 
not expect significant contributions from high redshift gas, for the 
purpose of our work we considered only the 25 snapshots in the 
interval < 2: < 1.5. 

The effect of a spatially unifor m, redshift dependent U V back- 
ground due to quasar radiation as in lHaardt & Madaul ( Il996 ) is con- 
sidered (with helium heating rates enhanced by a factor of 3, see the 
discussion in lTomatore et alj201(ll) . This process is important when 
considering the evolution of the amount of matter in the WHIM 
phase (see, e.g., Dave et al. 20()l l). Star formation is m o deled fol- 
lowing the multiphase model of ISpringel & Hemguistl ( |2003[) : in 
this prescription the SPH particles exceeding the density of nH=0. 1 
cm~3 are considered to host a hot phase and a cold one in pres- 
sure equilibrium, with the latter assumed to host the star-forming 
regions. Every two-phase SPH particle spawns up to three genera- 
tion of collisionless star particles in a stochastic way, each of them 
having about one third of the mass of the original gas one. 

During this proce ss the SPH part i cles a re enriched follow- 
ing the prescription of iTomatore et al.l ( l2007h that considers the 
formation of SN-II, SN-Ia together with low and intermediate 
mass stars i n the AGB phase. T h e corr e sponding stellar yields are 
tak en from 'Woosley & Weaver ('1995"), ' Thielemann et alj j2003l) 
and lvan d en Hoek & Groenewegen ( 1997) for the SN-II, SN-Ia and 
AGB, respectively. This code allows to follow independently the 
evolution of H, He and of 6 metal species, namely C, O, Mg, Si, 
S and Fe for every SPH particle. The formation of all the other 
elements is also computed in every gas particle though only glob- 
ally, storing their total mass in a single variable in order to save 
memory and computatio n time. The IMF adopted is the one pro- 
posed bv iKroupal j200lh . with a multi-slope function of the form 
(p(m) oc m~^, with the index y varying in the range 0.3-1.7 from 
low to high masses. 



Radiative cooling due to free-free interaction of the electrons 
and ions is computed using a primordial mix of H and He. In ad- 
dition to that, coolin g from heavy elements is followed adopting 
the cooling rates of S utherland &DoDitalll993h : in this scheme the 
global metallicity of every SPH particle, and not the individual el- 
ements, is considered assuming solar yields. However, as we will 
explain in detail in Section [T2l we stress that for what concerns the 
results presented in this paper our mock emission s pectra are com- 
puted considering all the elements provided by the lTomatore et all 
( 2007) code separately. We point out also that this cooling scheme 
assumes implicitly that the IGM is in collisional ionization equilib- 
rium: although this is not completely true, this can be considered 
as a safe approximation of the gas ionization state in the densest 
regions of the WHIM that dominate the emission. 

Galactic out flows are modeled following 

ISpringel & Hemguistl ( |2003|) . Winds are assumed to blow from 
the star-forming regions with a mass loading rate M„ = rjM^,, 
where A'h is the star-formation rate of each multi-phase particle 
and the efficiency is a free parameter (set to 2). SPH star-forming 
particles are then stochastically selected to become wind particles, 
with a velocity of ?;w=500 km/s and are decoupled from the others 
allowing them to travel freely for 10 kpc. This allows the wind 
particles not to get stuck inside the dense star-forming regions. 

In order to test the impact of the different parameters and as- 
sumptions of o ur reference W model, w e analyze also the results of 
other 7 runs o f iTomatore etal] ilOUt) . It is important to note that 
all the simulations were run assuming the same initial conditions 
so that they actually follow the evolution of the same volume and 
matter (with the exception of W75,5i2 that follows a larger volume) 
thus allowing a direct comparison between each other. In the fol- 
lowing subsections we will describe these runs focusing on their 
differences with respect to the reference model. 

2.1 Simulation parameters 

The W75.512 and W37,4oo runs assume the same reference model 
but with different volume size and resolution of the particles. These 
simulations can be used to assess how the limited volume and finite 
mass resolution affect our results. 

The W75.512 run has a box size of L =75 comoving 
Mpc. The number of particles is augmented to 2x5123 in order 
to maintain the same resolution. On the other side the W37,4oo 
run follows the same volume but with higher resolution. The num- 
ber of particles is increased to 2x400'', achieving a resolution of 
m-dm = 5.0 X W'^h-^ Mq and mgas = 9.4 x lO'' ft"^ Mq. The 
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softening is also reduced to e = 4.8 h^^ comoving kpc at 2; > 2 
and e — 1.6 physical kpc at z < 2. 

The number of snapshots available for these runs is reduced. 
In the interval < z < 1.5 we have 6 snapshots for W75,5i2 
and 16 for W37,4oo. However, our scheme for the reconstruction of 
the observed light-cone (see Section l3Tt allows us to minimize the 
impact of these differences. 



2.2 Initial mass functions 

The Way and Ws runs differ from the r eference model as for the 
IM F adopted , respe ctively the one from lArimoto & Yoshiil ( Il987h 
and lSalpeted ( Il955h . Changing the IMF implies a modification of 
the proportion of the SN types which reflects into different element 
ratios: top-heavier models produce a relatively higher number of 
SN-II and therefore are expected to provide more oxygen with re- 
spect to iron and carbon. 

Indeed, the Way run shows a global value of [C/O]=-0.3 for 
the WHIM at 2 = while the W and Ws runs have values close 
to -0.15. All the three models show a similar behavior of the red- 
shift evolution with [C/0] ratios growing with time due to the in- 
creasing contribution of SN-Ia (see the corresponding discussion in 
iTomatore et"ai]|201(]h . 



2.3 Feedback models 

The NW, CW and BH runs adopt a modified feedback scheme with 
respect to our reference W run. In the NW run galactic winds are 
simply turned off, while the CW simulation keeps wind particles 
hydrodynamically coupled with the other ones (more details about 
th e effect of coupled winds on s tar-forming systems can be found 
in lDalla Vecchia & Schav3l2008b . 

The BH run assumes no winds but includes the energy feed- 
back associated to supermassive black-hol es (SMBH) following 
ISpringell ( |2005|) and lDi Matteo et al.l ( l2005h . In brief, SMBHs are 
computed as sink particles hosted by haloes of mass higher than 
Afth = 10^" h~^MQ that grow with a rate of 

Mbh = min(MB, MEdd) , (1) 

where Mb an d Msdd are the Bondi-Hoyle-Lyttleton accretion 
rate (see, e.g.. lBondilll952 ') and the Eddington rate, respectively. 
The energy rate injected into the surrounding gas is computed as 
^jfccd = ei-efAfljji being and ef the two free parameters repre- 
senting the radiative efficiency and the fraction of energy coupled to 
the gas, res pectively. These values are set to = 0.1 and ef = 0.05 
with which loi Matteo et al.l ( l2005h were able to reproduce the ob- 
served relation between BH mass and velocity dispersion of the 
bulge. 

We refer the reader to the detailed discussion in 
iTomatore et al.l feoiof) about the effects of these models on 
the properties of the IGM. As a general remark, here we point 
out that NW and BH models show a stronger star formation 
rate (SFR) at higher redshift with respect to our reference run. 
Starting at z ~ 3 the BH feedback causes a heavy suppression of 
star-formation, getting closer to the W simulations, while for the 
NW the SFR stays higher even at later epochs. It is important to 
highlight also that BH feedback is able to affect significantly the 
properties of collapsed objects causing the IGM to be more diffuse 
at all scales. 



3 THE SIMULATED LIGHT-CONES 

Extracting observables from the outputs of LSS simulations re- 
quires a non-trivial post-processing analysis: indeed, one needs to 
reconstruct all the volume of the past light-cone as seen by an ob- 
server located at z = in order to account for all the matter that 
contributes to the integrated signal up to a given redshift. 

For the purpose of estimating the properties of the WHIM 
emission, we decide to consider the interval < z < 1.5: 
this is eno ugh to account for most of the expected signal (see 
iRoncarellTet al. 2006; Takei et al. 201 1). Then we compute the X- 
ray emission spectra associated to the gas elements present in the 
light-cone volume, we integrate them along the line-of-sight and, 
finally, we produce mock 3D-maps in which every pixel contains 
the spectrum as seen by a telescope. 

In this Section we describe in detail the different steps of this 
procedure. 



3.1 Reconstructing the light-cone volume 

The method desc ribed here is an improve ment of the one adopted in 
iRoncarelli et al.l 82006) and iUrsino et al.l ( I201Q) for the X-ray emis- 
sion of t he LSS (see also similar applications for estimates of the SZ 
effect in 'da Silva et al."200(l I2OO1I : IWhite et alj|2003: IZhang et al.1 
[2002; Roncarelli et al. 2007]), since it includes a more general ran- 
domization scheme and a more preci se procedure for the redshift 
sampling (as in , Roncarelli et al.ll2010h . 

We reconstruct the light-cone volume by replicating the out- 
puts of the simulations along the line of sight up to the redshift 
limit of z = 1.5. With the cosmological parameters assumed 
here (see Section [2} this corresponds to a comoving distance of 
rfc.iim = 3206 /i~^Mpc, which is 86 times the length of our simu- 
lation boxes (43 for the W75,5i2 run). When considering large dis- 
tances one needs also to take into account the evolution of structure 
formation by using the different simulation snapshots. To this pur- 
pose for a given distance from the observer we chose the output that 
better approximates the age of the Universe at that epoch. Since in 
our case we have more cubes to stack than available outputs, most 
of the cubes will correspond to the same snapshot than their neigh- 
bor(s). When a snapshot change intervenes the cube is split into two 
slices at the corresponding distance. 

In order to avoid the repetition of the same structures in our 
maps, every cube undergoes a series of randomization procedures. 
First we give to every axis a 50 per cent probability of being re- 
flected, then we chose randomly the axis order. Once this is done, 
we exploit the periodic boundary conditions of our simulations to 
select a position inside the cube which is then put in its center. Fi- 
nally we introduce a rotation of a random angle along the line of 
sight only if the distance from the observer is small enough to en- 
sure that the projected cube encloses all the map surface. As in our 
previous works, the slices corresponding to different snapshots but 
belonging to the same cube undergo the same randomization pro- 
cess to preserve power on large scales. 



3.2 Computing the X-ray spectra and maps 

For every SPH particle that falls inside the light-cone we compute 
its spectrum deriving it from its physical variables. We adopt the 
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APEC emission model jSmith et al.ll200ll . version 2.(0) which as- 
sumes that the emitting gas is an optically thin plasma in colli- 
sional ionization equilibrium. We create a series of emission tables 
by using the VAPEC implementation on XSPEC (version 12.6.0) that 
allows to separate the contribution of the different elements. More 
in detail, for every element considered in our simulations (H, He, 
C, O, Mg, Si, S, Fe and the sum of all the other|fl as described 
in Section|2]( we extract three-dimensional tables Sci(z, T, E) with 
fixed metallicity and normalization in which two variables are the 
redshift and temperature and the third one is the energy channel. 
Then the flux as a function of energy E is computed for the i-th 
SPH particle as 

F^iE) ^ ""^'l^y^' ^Xei,s.i(^.,r.,j;) , (2) 

where nc,i and nu.i are the electron and hydrogen number density, 
respectively, for each particle, Vi is its physical volume and dc.i is 
its comoving distance from the observer. The redshift Zi and the 
temperature Ti of the SPH particle are used to interpolate the value 
of Sci{zi, Ti, E) from the tables. The abundance of each element is 
defined as 

= Aci — '- , (3) 

where Ad is the element atomic mass in units of the proton mass 
and nci,i is the number density of the element provided by the 
codell 

Finally, in order to obtain the surface brightness the flux of 
every particle computed from eq. ^ is distributed in the map pixels 
by integrating analytically the SPH smoothing ker nel (the details of 
the procedure are explained in lUrsino et alj2010h . 



which we will make an extensive comparison. In practice, the dif- 
ference consists in the inclusion or not of the warm-hot dense gas 
that roughly corresponds to galaxy groups. As a reference, with our 
definition the WHIM accounts for 33 per cent of the total baryonic 
mass at 2 = in the W run. 

Given the redshift limit of 2: < 1.5, the maximum angular 
extension of the maps is fixed by the box length and corresponds 
to about 40 (80 for the W75,5i2 run) arcmin, however we chose to 
compute maps of 30 arcmin per side in order to avoid problems 
caused by border effects: this size roughly corresponds to the dou- 
ble of the CDFs, so that every mock map covers four times as much 
sky surface. 

The first two sets of mock maps, as indicated in Table |2] have 
512 pixels per side which corresponds to a resolution of 3.52 arc- 
sec. The spectra span the energy range 0.3-10 keV and are sampled 
with bins of 50 eV for energies lower than 2 keV. At higher ener- 
gies we only compute the surface brightness of four 2 keV-wide 
bands. The results of these two sets will be compared mainly with 
Chandra and XMM-Newton data in the following section. 

Our method of spectra computation, as described by eq. l|2]l, 
allows us to separate the contribution of the different chemical ele- 
ments. Taking advantage of this, we can artificially "switch off" the 
emission of all other elements in order to focus on the properties of 
oxygen only. To this purpose we compute two other sets of mock 
maps that include only the signal from this element. Since our goal 
is to study line emission, we increase the spectral resolution to 1 
eV, while reducing the energy range to an upper limit of 1 keV and 
the angular resolution to ~ 1 arcmin. We produce other two sets 
(for the IGM and the WHIM) with this method and use it to predict 
oxygen line counts (see Section|5j- 



3.3 The light-cone set 

The method described above allows us to compute a set of 3-d mock 
spectral maps (i.e. the surface brightness SB{0,E), as a function 
of the position in the sky 9 and energy E) from which we extract 
our results. 

To assess the statistical robustness of our results we construct 
20 different light-cones by varying the initial numerical seeds nec- 
essary for the randomization process (see Section 13.11 ). Each of 
them is applied to all of the simulation runs, ensuring that every 
light-cone represents the same comoving volume (except for the 
W75,5i2 run), simulated with different physical prescriptions. This 
is important because it allows us to compare directly the results 
of our simulations without dealing with uncertainties associated to 
cosmic variance. On the whole we obtain 640 outputs with different 
characteristics summarized in Table [2] 

In our analysis we distinguish between the properties of the 
IGM as a whole and of the WHIM. Here we adopt the definition 
of WHIM as the gas with temperature in the range 10^-10*^ K and 
density contrast S < 1000. We advise the reader that not all the 
authors adopt the same definition: in some works the WHIM is de- 
fined onl y with respect to the te mperat ure cut, in particular th is ap- 
plies to jRoncarelli et al.) ^20061) and to lTomatore et al.l ( I2OI0I) with 

^ We are aware of the error in the calculation of some ions present in this 
version of the code. We checked that its impact on our results is negligible. 
^ The relative a bundances of the ele ments not considered individually are 
computed from Asplund et alj i2009l) . 

^ In the formula of eq. J2) the abundance is generalized also for hydrogen 
for which, of course, Xjj = 1- 



4 PROPERTIES OF THE X-RAY EMISSION 

We show in Fig. [T] the expected IGM surface brightness in the 
0.3-0.8 keV band for a single light-cone with 4 different feedback 
schemes. As said in the previous section, with our volume recon- 
struction procedure we ensure that our maps reproduce exactly the 
same structures (i.e. the same objects in the same positions), with 
no difference associated to cosmic variance. 

The emission is dominated by galaxy clusters, with peaks that 
have similar values for all the models (~1 ph s~^ cm~^ deg~^). 
A significant fraction of the maps has a fainter surface brightness 
from 10~^ down to 10~^ ph s~^ cm~^ deg~^. In the W model the 
haloes are more diffuse due to galactic winds that eject gas from 
the centre: given the relatively small size of our simulation vol- 
ume we should consider this as representative of the effect at the 
galaxy group scale. When turning winds off (NW model) clusters 
have sharper peaks, while the CW model produces an intermediate 
effect. In the BH model the cluster emission is significantly lower: 
the strong energy injection at early epochs increases considerably 
the entropy of the gas causing a suppression of gas accretion from 
high redshift. As a result smaller haloes are disrupted while for big- 
ger ones only the ce ntral peaks are present (see the discussion in 
iTornatore et al.l2o"lol) . On the contrary, with this model low surface 
brightness regions are slightly brighter: apart from cluster centers 
all pixels stay roughly within 2-3 orders of magnitudes in surface 
brightness. Due to the very strong effect on galaxy clusters and 
groups, this model has to be considered as an extreme and non- 
realistic case. It is however useful as a case of study to understand 
the effect of this kind of feedback. 

The maps of Fig. |2] show the surface brightness of the same 
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Figure 1. Surface brightness of the IGM (i.e. all the gas) in the 0.3-0.8 keV band. The maps are 0.5 deg per side and represent the same light cone simulated 
with different feedback schemes: galactic winds (W model, top left), no winds (NW model, top right), coupled winds (CW model, bottom left) and black-holes 
feedback (BH model, bottom right). 
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Figure 2. The same as Fig.[T]but for the WHIM (i.e. gas with lO'^ K< T < lO'' K and (5 < 1000). 
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Table 2. Summary of the sets of 3-d spectral maps created for our studies. In all cases the redshift interval is < 2 < 1.5. Column 1: total number of mock 
maps (the number of simulation runs times the number of independent light-cones). Column 2: gas phase, IGM corresponds to all of the gas present in the 
simulation outputs, WHIM con'esponds to the gas with temperature in the range 10" — lO*^ K and density contrast 5 < 1000. Column 3: chemical elements 
considered. Column 4: angular size if the maps. Column 5: angular resolution. Column 6: energy range of the spectrum. Column 7: energy resolution for the 
E <1 keV bins. 



TV Gas phase Elements Map size Ang. resolution Energy range En. res. (E < \ keV) 



(sims X cones) 






(deg) 


(arcsec) 


(keV) 


(eV) 


8 X 20 


IGM (all) 


All 


0.5 


3.52 


0.3-10 


50 


8 X 20 


WHIM 


All 


0.5 


3.52 


0.3-10 


50 


8 X 20 


IGM (all) 


O 


0.5 


56.25 


0.3-1 


1 


8 X 20 


WHIM 


O 


0.5 


56.25 


0.3-1 


1 



light-cones but considering the WHIM only: the difference with re- 
spect to the previous ones is of course the absence of the hot bright 
cluster centers. These images show that even if our definition of 
WHIM includes the density cut at <5 < 1000, cluster and group 
outskirts close to the virialized regions are still present. However, 
these structures can be considered below the detection limit of the 
CDFs (see Section|4T|and Fig.lU. It is important to note that even 
when excluding artificially the hot gas directly from the simula- 
tions, no filamentary structure is clearly visible in these maps. This 
is due to the superimposition of the structures of the cosmic web 
along the line of sight with the effect of self-washing out. Also in 
these maps the most important differences are in the BH model, 
where the WHIM appears more diffuse and fainter than in the other 
cases. The effect of enhancing the surface brightness of the faintest 
pixels has little impact on the possibility of detecting these regions 
since they still fall well below the limits of current instruments. 
The images of the other models (not shown in these figures) do not 
differ significantly from our reference one (W model). 

We analyse the effect of feedback on the spectral properties of 
the IGM and the WHIM with the plots shown in Fig. [3] and with 
the corresponding surface brightness in different bands shown in 
Table |3] When considering a large box (W 75, 512, top panel), the 
global IGM emission is about 20 per cent higher at all energies 
due to the presence of more massive clusters. On the contrary, no 
significant change is produced when considering the WHIM emis- 
sion, as an indication that our box size is large enough to sample 
the LSS perturbations that influence the WHIM physics. The op- 
posite effect applies when increasing the resolution of our simula- 
tions (W37,4oo): almost no difference can be seen for the IGM sig- 
nal while the WHIM emission increases by ~10 per cent. This is 
mainly due to the effect of resolution on metal formation/emission, 
as we will discuss in Section |5] On the whole we can safely con- 
clude that the limited box size and resolution of our simulations 
have a negligible impact on the expected surface brightness of the 
WHIM. 

When considering different IMFs (Way, Ws, central panel), 
the most interesting effect to point out is the fact that the Way 
run produces a lower surface brightness for the IGM, dominated 
by T > 10^ K gas, and a higher one for the WHIM. Since the 
lArimoto & Yoshii ( [l987l) IMF produces a higher amount of mas- 
sive stars, this turns into a relatively lower amount of Fe, produced 
by SN-Ia, and a correspondingly higher amount of O eje c ted by 
SN-II. The opposite applies when assuming the ISalpeteJ jl955l) 
IMF. Since O has the peak of its emission lines at T ~ 10® K 
and energies around 0.5 keV, the maximum of the difference can 
be seen for the surface brightness in the softest 0.3-0.8 keV band. 



in which the WHIM is 20 per cent brighter for the Way run with 
respect to the Ws one. 

The most important variations are introduced when consid- 
ering the different feedback schemes (bottom panel). On one side 
galactic winds suppress the peak emission in the central regions of 
clusters, while at the same time they allow the gas to escape the 
cores and to produce significant emission also in the outskirts (see 
Fig.lTJ. As a result the IGM surface brightness is increased by a fac- 
tor of ~ 2 at all energies when comparing the NW and W runs, with 
the CW run that constitutes an intermediate scenario. The same ap- 
plies to the WHIM emission but with only a ~ 20 per cent increase 
from the NW to the W run. This is particularly evident at i5 < 0.8 
keV where oxygen emission dominates, indicating that winds are 
able to spread metals in the regions at intermediate temperatures. 

On the other side, BH feedback greatly suppresses the emis- 
sion, for both IGM and WHIM. In the IGM the difference range 
from a factor of 5 at lower energies up to one order of magnitude in 
the 2-10 keV band, while for the WHIM it is a factor of ~ 3. This is 
the result of the strong BH winds that prevent gas from collapsing 
and suppress star formation. 

We also estimate the cosmic variance associated to our results 
by computing the standard deviation (quoted as error in Table[3]( in 
80 fields of 225 arcmin^ . For the IGM the cosmic variance is very 
high, comparable to the value itself: this happens because the pres- 
ence of clusters and groups in the field dominates the total signal. 
Since the WHIM is dominated by the diffuse component its rela- 
tive cosmic variance is lower, being of the order of 30 per cent in 
the 0.3-0.8 keV band, but still higher than the expected differences 
introduced by the various physical prescriptions. 

We can compare our results for the IGIvQ also with our pre- 
vious work tRoncarelli et al.u2006) . Even if the feedback model 
of our W run is s imilar to the one adopted in that work (see 
iBorgani et al.ll2004l) . their total expected emission in the 0.5-2 and 
2-10 keV bands is higher by a factor of 2 and 3, respectively. This 
difference is due to the larger volume of the simulation that allowed 
to better sample the high-mass tail of the mass distribution, thus ob- 
taining brighter haloes that dominate the emission. 

When c ompa ring our WHIM results with the ones of 
lUrsino et al. I ( l2010() in the 0.65-1 keV band, our predictions are 
close to their 0.15 x 10~^^ erg s^^ cm^^ deg^^ obtained with 
a semianalytical characterization of the metallicity and a factor of 



* We do not compare our results on the WHIM given the different defini- 
tions (with and without the density cut at 5 < 1000) adopted in the two 
works. 
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Table 3. Expected surface brightness (in units of 10 erg s ^ cm ^ deg ^) for the whole IGM (i.e. all the gas) and for the WHIM (i.e. gas with 
T = 10^ - 10"^ K and 5 < 1000) in different X -ray bands. In each line the results ai'e shown for the 8 different simulations as the mean computed over 20 
independent light-cones (a total of 5 deg^). The errors represent the standard deviations in fields equal to a quarter of the original maps (e.g. 225 arcmin^, 
close to the size of the CDFs). 









Surface brightness (10 er; 


— 1 —2 

5 s cm 


deg- 2) 








IGM 


WHIM 


IGM 


WHIM 


IGM 


WHIM 


IGM 


WHIM 


Model 


(0.3-0. J 


i)keV 


(0.65- 


-1) keV 


(0.5-2) keV 


(2-10) keV 


W 


1.05±1.08 


0.29±0.10 


0.73±1.14 


0.12±0.09 


2.12±3.84 


0.28±0.16 


0.52±1.84 


<0.01 


W75,512 


1.20±0.90 


0.30±0.13 


0.84±0.86 


0.13±0.10 


2.45±2.72 


0.30±0.18 


0.61±1.28 


<0.01 


W37,400 


1.07±1.18 


0.33±0.11 


0.76±1.30 


0.14±0.09 


2.18±4.29 


0.32±0.17 


0.52±1.94 


<0.01 


Way 


0.94±0.90 


0.31 ±0.11 


0.64±0.94 


0.13±0.10 


1.88±3.22 


0.30±0.17 


0.48±1.57 


<0.01 


Ws 


1.17±1.08 


0.26±0.09 


0.82±1.14 


0.10±0.07 


2.32±3.34 


0.24±0.13 


0.46±1.03 


<0.01 


NW 


0.57±0.50 


0.24±0.08 


0.38±0.53 


0.11 ±0.08 


1.06±1.84 


0.24±0.13 


0.26±0.94 


<0.01 


cw 


0.84±0.71 


0.28±0.10 


0.55±0.73 


0.12±0.09 


1.59±2.53 


0.28±0.16 


0.39±1.27 


<0.01 


BH 


0.22±0.12 


0.10±0.03 


0.11±0.13 


0.04±0.03 


0.31 ±0.40 


0.09±0.05 


0.06±0.16 


<0.01 



^ 2 lower than their higher metallicity models. Even if these re- 
sults can be considered to be in a broad agreement, their differ- 
ences reflect the current uncertainties in how metals spread in the 
lower density regions and on the ability of hydrodynamical models 
to provide precise predictions on the WHIM emission. More pre- 
cisely, the typical metal abundances associated to our WHIM phase 
are of the order of Z ^ 0.05 Z©, that corresponds to a regime 
where the expected signal is roughly equally distributed between 
free-free and line-emission. When assuming a higher metallicity 
(Z — 0.1 — 0.15 Zq) the emission becomes metal-dominated and 
increases significantly, thus justifyin g the differences in t he results 
with the higher metallicity models of lUrsino et al.l bOld) . 



4.1 Comparison with the unresolved X-ray background 

Our results can be compared with the upper limits provided by the 
unreso lved X-ray background (UXRB) from lHickox & MarkevitchI 
( |2007|) . who measured the diffuse emission in the CDFs after sub- 
tracting all possible contaminating sources. However, we must con- 
sider that a direct comparison between our IGM maps/spectra and 
the measured UXRB is not possible since our mock maps include 
also clusters and groups that have been identified and eliminated in 
the CDFs observations. On the other side, using our WHIM maps 
can lead to misleading results since our WHIM definition is based 
on the intrinsic thermodynamical properties of the SPH particles, 
that do not correspond to an observ ational criterium of sou rce iden- 
tification (see also the discussion in lRoncarelli et al.ll2006l ). 

For these reasons we decide to adopt an observationally ori- 
ented approach to estimate the UXRB that corresponds to our sim- 
ulations. Starting from our IGM maps/spectra (e.g., Fig.O we com- 
pute the photon counts in the 0.5-2 keV band expected from a 
Chandra Deep Field South (CDF-S) observation, by considering 
the ACIS-I response function at different energies (sampled in bins 
of 50 eV each). In order to simulate reliable Chandra images we 
added an isotropic level of background that takes into account the 
contribution of both galactic and particle background. We neglected 
the contribution of undetected AGN/galaxies since, at flux limits of 
the CDF-S, their flux can be considered completely resolved into 
point sources. The images were then corrected for an average ex- 
posure of 4 Ms with an artificially added Poisson noise. The source 
identification has been performed by running a simple sliding cell 
detection with a cell size of 12 pixels and fixing a signal-to-noise 



ratio threshold of 4. We have then excluded all the regions within 
which the overall encircled signal from sources is above 4a with 
respect to the background. An example of the implementation of 
this method is shown in Fig.|4l 

Once we have identified all the pixels corresponding to de- 
tected haloes, we compute the expected spectrum of the remaining 
regions for all the different light-cones and simulations. We show 
in Fig. |5] our results (averaged over the 20 light cones) for the W 
model for both the IGM a nd the WHIM and compare the m with the 
observational estimates of lHickox & Markevitc H j2007h . The same 
comparison is shown in Table |4] for the emission in three different 
X-ray bands and for the other physical prescriptions. 

When looking at the results for our reference model, the sur- 
face brightness per unit energy of the UXRB associated to the IGM 
is well represented by a power-law (SB« E~~') with spectral index 
7 =0.9 (1.5, for the WHIM) in the 0.3-0.8 keV band, with a clear 
steepening at higher energies, up to 7 =2.5 (3.8) in the 0.8-2 keV 
band. These values do not change significantly when other models 
are considered. 

Through all the energy range our IGM predictions fall below 
the upper limits of the UXRB. At energies E < 0.7 keV the ex- 
pected emission of the whole IGM is more than an order of magni- 
tude below the observed UXRB. This indicates that at soft energies 
other sources are the main contributors, likely the emission asso- 
ciated to local components (Galaxy and Solar System). For these 
reasons the search for LSS emission in the UXRB at these ener- 
gies requires a subtraction of these foreground components with 
more than 10 per cent accuracy. A possible observational technique 
consists in relying on their different angular correlation of the two 
components, as we will discuss in Section |4~2l 

On the contrary, we predict that at energies E > 1 keV a 
significant fraction (from 50 per cent to all) of the UXRB is due 
to the emission of the IGM. The high value of the cosmic variance 
(reported as error in Table|4l( together with the low angular surface 
of the CDFs does not allow us to provide more precise estimates. 

When comparing our models with the different feedback and 
IMFs assumptions, the same considerations apply with respect to 
the discussion on Table [3] galactic winds increase consistently the 
expected IGM signal. 

It is important to point out that in all cases the surface bright- 
ness associated to a diffuse WHIM component is less than half of 
the total IGM one, ranging from --^20 per cent to --^40 per cent 
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Table 4. The same as Table [3] but after subtracting clusters and groups identified in a mock CDF-S observation. The last line indicates the observational 
measurements of .Hicko x & Markevitch (2007) with relative errorbars. 







Surface brightness (10 


12 — 1 - 

erg s cm 


2 deg-2) 






IGM 


WHIM 


IGM 


WHIM 


IGM 


WHIM 


Model 


(0.5-2) keV 


(0.65- 


-l)keV 


(1-2) keV 


W 


0.57±1.31 


0.17±0.21 


0.23±0.54 


0.08±0.12 


0.23±0.67 


0.04±0.06 


Way 


0.45±0.77 


0.19±0.24 


0.18±0.34 


0.09±0.13 


0.17±0.37 


0.04±0.06 


Ws 


0.48±1.05 


0.14±0.18 


0.20±0.49 


0.06±0.09 


0.18±0.46 


0.04±0.06 


NW 


0.33±0.42 


0.16±0.17 


0.14±0.19 


0.08±0.10 


0.11 ±0.20 


0.03±0.04 


CW 


0.41±0.69 


0.17±0.22 


0.17±0.31 


0.08±0.13 


0.14±0.32 


0.04±0.06 


BH 


0.20±0.18 


0.08±0.05 


0.08±0.07 


0.04±0.03 


0.07±0.09 


0.02±0.01 


Obs. 


4.3±0.7 


1.0±0.2 


0.34±0.14 



from high to low energies, with the only exception of the 50 per 
cent of the Way model in the 0.65-1 keV band. This suggests that 
a consistent part of this signal is not associated with truly diffuse 
(5 < 1000) gas. In particular this becomes important at i5 > 1 
keV where the majority of the UXRB due to the LSS is clearly 
associated to cluster outskirts and faint galaxy groups. However 
we must consider that LSS filaments also contain dense clumps 
that are excluded by our WHIM definition so, as a bottom line, we 
can consider our two estimates, IGM and WHIM, as bracketing the 
emission associated to the missing baryons component. 

4.2 The angular correlation function 

Given the relatively low surface brightness of the WHIM, one of the 
possible techniques to disentangle its signal from the other brighter 
components (e.g. Galactic emission) consists in performing a cor- 
relation function analysis. To this purpose we compute the angular 
(auto)correlation function {AcF) w{6) of the surface brightness in 
the 0.3-0.8 keV band. 

We proceed as follows. First, for a given map we compute its 
surface brightness contrast 

5{x) = {SUx)/S^) - 1 , (4) 

where Sh{x) is the surface brightness at the position x and Sh is 
the average of the map. Then we compute the angular correlation 
function as 

w{e) ^ {S{x)5{x + 6)) . (5) 

With this definition a value of w{6) = corresponds to a random 
uncorrelated field. 

We show in Fig. [6] our results for the IGM and the WHIM 
and for our 4 different feedback scheme s. As expected (see, e.g., 
iRoncarelli et al Jl 20061 ; [Ursino et al.]|201 ih . the IGM correlates sig- 
nificantly up to angular scales of 2-3 arcmin, which correspond to 
the typical size of galaxy clusters and groups. At larger scales it 
drops rapidly to values consistent with zero (9 > 5 arcmin). Again, 
the cosmic variance produces significant deviations (grey-shaded 
area) from map to map. Galactic winds enhance the correlation sig- 
nal of the IGM due to higher relative clusters emission. The BH 
model produces a high correlation only at very low angular scales 
(9 < 0.2 arcmin), consistently with what appears from Fig.[T] The 
WHIM AcF is lower by almost one order of magnitude at low an- 
gular scales but has a shallower slope making it comparable to the 
IGM one for 9^2 arcmin. On the other side, the different wind im- 



plementations (W, NW, CW) produce small variations in the AcF, 

of the order of 10 per cent. 

We can compare our results with lGaleazzi et al] | |2009|) who 
detected a putative WHIM emission by analysing the AcF of 6 
different XMM-Newton fields in the 0.4-0.6 keV bancfl In this 
comparison one must take into account the different normalization 
adopted with respect to our work: in fact, since they do not know 
a priori how much emission is associated to a diffuse component, 
they normalize w{6) to the total UXRB flux, including also the 
foreground contribution. We can obtain an order-of-magnitude es- 
timate of the impact of this difference by using our results on the 
average 0.4-0.6 keV surface brightness associated to unresolved 
IGM emission (see Section l 4. U an d compare it to their total UXRB 
from lHickox & Markevitc ij j2007h : the ratio between the latter and 
the former is ~ 26, which translates into a difference in the nor- 
malization of the order of ~ 70 0. If we apply this correction to 
the data of iGaleazzi elal] j2009h . we obtain a value of 'w{9) ~ 7 
for 9=\-2 arcmin, that falls within a factor of ~ 3 from our IGM 
results. Given the high values of relative errors and cosmic vari- 
ance associated to these measurements, we can not obtain a more 
precise estimate, but we can c onclude that o u r resu lts appear in 
broad agreement with those by IGaleazzi et al.l (|2009|) at small an- 
gular scales. On the other side, we obtain a steeper slope of the 
AcF at large angular scales, suggesting the possible presence of a 
contribution from other components in their data. 



5 OXYGEN LINE EMISSION 

In this Section we focus on the properties of oxygen emission. In 
particular, we study the expected surface brightness due to Ovil 
and OviII that represent the ions with the most significant contri- 
bution at X-ray energies. 

To this purpose we take advantage of our method for comput- 
ing the spectra associated to the SPH particles described by eq. ^ 
that allows us to separate the emission from the different elements. 
As discussed in Section [331 we use the last two light-cones set 
listed in Table |2l which were computed with an energy resolution 
of 1 eV necessary for the analysis of single line emission. The en- 
ergy range considered here (0.3-1 keV) allows us to enclose all of 



^ We checked that the changes in the AcF associated to the two different 
bands are negligible. 
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Figure 3. Expected spectra of the IGM (solid lines) and WHIM (dashed 
lines), computed as average over the 20 independent light-cones. In all 
plots the blue lines correspond to our reference model (W). Top panel: com- 
parison with W75 512 (green) and W37 400 (red). Central panel: compar- 
ison with Way (green) and Ws (red). Bottom panel: comparison with NW 
(green), CW (orange) and BH (red). 




Figure 4. Sketch of the haloes (marked by the elliptical regions) identified 
in one of our mock maps. The map represents the surface brightness (on a 
logarithmic scale) at 0.3 keV and corresponds to the same light cone shown 
in the maps ofFigs.Handll] 
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Figure 5. Expected spectrum of the diffuse XRB associated to extragalactic 
gas emission for the W model. The solid line corresponds to the whole 
(IGM) gas emission, the dashed line to the WHIM contribution. The lines 
correspond to the median relative to 80 different values computed on an 
area of 225 arcmin^ (e.g. CDFs), while the shaded region encloses the first 
and third quartile for t he IG M. The black squares represent the data from 
iHickox & MarkevitchI i2007l) with corresponding errorbars. 



the most important OVII and OviII lines up to z ~ 0.9, with no sig- 
nificant contribution from Ovi. In order to save memory, these O 
emission spectra have been computed with a poorer angular resolu- 
tion of ~ 1 arcmin, but this does not introduce any systematics since 
the expected correlation of the IGM and WHIM signal does not rise 
significantly at lower angular scales (see Section l4~2t . Moreover the 
expected angular resolution of the next-generation spectrographs 
(e.g., the IXO/ATHENA proposaQ) is higher than this limit. An ex- 



|http://sci.esa.int/science-e/www/object/index.cfm?fobjectid=48729| 



ample of a WHIM spectrum computed with our method for both 
oxygen-line and continuum emission is shown in Fig. [7] 

We compute the statistics of the expected number of O lines 
(OviI and OVIII added together) above a given surface brightness 
along the line of sight. We proceed in the following way. For a 
given pixel of our maps, that corresponds to a line-of-sight of 0.88 
arcmin^ and extending in the range Q < z < 1.5, we consider its 
spectrum in the energy interval 0.3-1 keV and identify the 1 eV 
bins corresponding to a surface brightness higher than a threshold 
fixed to SBth = 2 X 10~^ ph s~^ cm~^ deg~^ (safely below the 
detection limits of future instruments). Then, since an emission line 
can extend to more than 1 eV, we assign to a single line all the 
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Figure 6. Angular correlation function of the IGM (solid lines) and the 
WHIM (dashed) for different feedback models: the reference W model 
(blue line), the NW model (green), the CW model (orange) and the BH 
one (red). Each curve represents the median computed over the 20 indepen- 
dent light-cones (area of 900 arcmin^ ) while the shaded regions enclose the 
first and third quartile. 
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Figure 7. Example of a spectrum extracted from a pixel of 1 arcmin^ size, 
for the WHIM component only. The black line shows the continuum emis- 
sion (i.e. including all the elements except oxygen) computed with 50 eV 
energy resolution. The green line represents the spectrum with the addition 
of oxygen emission, computed with 1 eV energy resolution to allow the 
identification of the lines. 

contiguous bins that fall above this threshold and sum their values 
to obtain the total surface brightness of the line itself. We repeat the 
same procedure for all the pixels and all the 20 different light-cones 
for both the IGM and the WHIM. We checked that the final results 
do not change significantly when degrading our angular resolution 
to ~2 arcmin or when changing the value of SBth. 

Fig. [8] shows the expected number counts of oxygen lines for 
our reference run compared with the other ones. When considering 
the different simulation parameters (top panel) we see that increas- 
ing the mass resolution of the simulation by a factor of ^^4 enhances 
the expected line counts for the IGM by ~20 per cent. This is due 
to improved capability of resolving the high density regions that 
host the star-forming metal-rich gas and are therefore bright line- 
emitters. The effect is also present when considering the WHIM 
oxygen lines in the faint (SB < 10~* ph s^^ cm^^ deg^^) end. 
This suggests that a high resolution is mandatory in this kind of 
analysis and is more desirable than a large simulation volume. 

The different IMFs (central panel) produce changes in the 
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Figure 8. Expected counts of oxygen (Ovil and OviII) line emission as- 
sociated to the IGM (solid lines) and WHIM (dashed lines) as a function 
of lines brightness. In all plots the blue lines correspond to our reference 
model (W). Top panel: comparison with W75,5i2 (green) and W37,4oo 
(red). Central panel: comparison with Way (green) and Ws (red). Bottom 
panel: comparison with NW (green), CW (orange) and BH (red). 

amount of SN-II and, consequently, on the total amount of oxy- 
gen released by stars. With the top-heaw lArimoto & Yoshiil jl987h 
IMF (Way model) the increased amount of SN-II produces signifi- 
cantly more oxygen. If we take as a reference a surface brightness 
detection threshold of 3xlO~^ ph s~^ cm~^ deg~^ (a possible 
goal of next-generation instruments), this translates into ~20 per 
cent more detected lines for both IGM and WH IM. The opposite 
effect happens when adopting the lSalpete3 ( Il955h IMF (Ws). 

The changes introduced by different feedback schemes (bot- 
tom panel) are also relevant. When suppressing galactic winds 
(NW) metals remain more concentrated in the cold star-forming 
regions where they are not able to emit significantly: this results in 
a factor of 2 less expected lines for the IGM. The CW model shows 
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trend for the IGM similar to thie W model, with a small increase for 
fainter lines. When looking at the WHIM curves, the effect of metal 
transport due to galactic winds can be seen with the higher amount 
of lines expected with respect to the NW model, that reaches an 
order of magnitude for brighter (SB> 10"* ph s~^ cm~'^ deg^^) 
lines. Coupled winds are not able to diffuse efficiently metals in 
the low density regions, resulting in a number of lines closer to 
the NW run than to the W one. Considering the BH feedback, the 
expected line counts are significantly smaller for both IGM and 
WHIM. Again, this is the result of the early energy feedback that 
suppresses the gas collapse resulting in the IGM having less dense 
structures. This effect is even more remarkable in the WHIM. 



5.1 Comparison with other models 

Apart from the differences introduced by the various feedback 
mechanisms, we can compare the results of our reference run with 
the expected oxygen li ne counts estimated by iTakeietalJ hoill) 
an d|c en & Fana ilOO &i. On the whole, the number of oxygen lines 
that we estimate is significantly lower with respec t to their results. 
If we consider the B 1 model of iTakei et al.l ( 1201 ih and a reference 



surface brightness of O.I ph s cm sr , they obtain a value 
of A'^ ~ 2 for the sum of Ovil and OviII lines in the IGM case 
and AT ~ 1 for the WHDvfl while our reference model predicts 
0.5 and 0.25 respectively, a factor 3^ lower. The reason of this 
discrepancy has to be searched in the di fferent recipe adop ted to 
describe the metal enrichment of the gas. iTakei et ^ j20I ll) adopt 
a post-processing method assigning a metallicity to the SPH par- 
ticles a s a function of overdensity: in particular, they use the re- 
sults o f lCen&Ostrikeilfr99a who estimated a relation of the form 
Z PC 6^^^. If we compare it to what shown in iTornatore et al.l 
1 I2OIO I) we can see that this results in a factor of 2 less metals at 
overdensities 5 < 10''. In addition to that, if we take into account 
also the differences in the relative eleme nt abundance, our model 
presents a 0/Fe mass ratio of ~ 2.5 while iTakei et all ( I20I 1*) adopt 
solar yields from I Anders & Grevessd ( Il989h which correspond to 
a mass ratio of 5.2. This adds another factor of 2 of differences in 
the abundance of oxygen. At a first approximation, given the lin- 
ear dependence of line emissivity on the element abundance, this 
translates into a factor of ~4 in line surface brightness that is suf- 
ficient to explain the differences with their work. In order to verify 
this hypothesis we run a test light-cone implementing their Z — S 
relatio n and, as expected, we obtained results similar to lTakei et al.l 

( I2OIII) . 

We can also compare our results with lCen & Fand ( |2006|) who 
computed the expected oxygen lines by using a hydrodynamical 
simulation. Even if they provide only the value of dN/dz at z = 0, 
we can est imate that our num bers are smaller, coherently to what 
said about Ir^kei et aLll l20Ill) . Again, part of the d ifference is due 
to the different element vields: ICen & Fan3 ( 120061) artificially aug- 
ment the oxygen abundance of their model (see the discussion in 
their Section 2) to reach a value of 0/Fe close to 2 times solaJf[ 
Such a high val ue is ruled out by current estimates from o bserved 
galaxy clusters dTamura et alj2004l : ISimionescu et alj2009h . In ad- 
dition to that, their simulation is able to spread more efficiently the 



^ The authors show the value of AN/Az. Since their redshift interval is 
< 2 < 0.5 here we consider that the contribution from lines associated 

to higher redshift emitters is negligible^ 

* We converted their 1.67 refer red to lAnders & Grevessd )l989l) yields to 
the ones of lAsplund eliai] j20O9l) 



metals outside the hot and dense regions, resulting in more oxygen 
in the thermodynamical conditions that correspond to bright OviI 
and OviII emission. This issue is connected to the current uncer- 
tainties in the feedback implementation and of the hydrodynamical 
codes themselves that require further investigations. In particular, 
an interesting comparison can be provided by estimates from Ovi 
absorption systems that, however, are beyond the scope of this pa- 
per. 



6 CONCLUSIONS 

In this paper we have analyz ed a set of cosmologica l hydrodynam- 
ical simulations presented in ITornatore et al.l ( I2OI0I) that include a 
detailed treatment of chemical enrichment of the IGM. The aim 
is to provide a description of the expected X-ray emission of the 
WHIM, with focus on the influence of different feedback mecha- 
nisms, namely initial mass functions, galactic winds and black hole 
feedback. 

Our results can be summarized as follows. 

(i) Galactic winds have a strong impact on the expected signal of 
the WHIM and the IGM in general: with our assumptions (v^ — 
500 km/s, 77 = 2) the surface brightness is enhanced by a factor of 

2 through all the energy range of interest (0.3-10 keV). This is due 
to the ability of winds to push gas and metals in cluster outskirts 
and outside, making lower density regions brighter. 

(ii) On the other side BHs, acting at earlier epochs, can prevent 
the IGM to collapse in dense structures suppressing significantly 
the expected emission of clusters and WHIM. 

(iii) All of our models predict that the expected contribution from 
non-resolved LSS structure s (galaxy groups, WHIM) is below the 
upper limits of the UXRB dHickox & Markevitc h' '20071 . At ener- 
gies E < 0.7 keV they account for only ~ 10 per cent of the 
UXRB, indicating that Galactic components dominate, while at 
E > 1 keV, they account from half to all of the signal, with signif- 
icant field-to-field differences expected due to cosmic variance. 

(iv) The WHIM shows an angular correlation function ■w{6) > 1 
up to scales of 2-3 arcmin, with very low variations due to feed- 
back effects or cosmic variance. At low angular scales (^ = 1 — 2 
ar cmin), our resu l ts are in broad agreement with the observations 
of iGaleazzi et al.l ( l2009h . while for larger angles we obtain lower 
values. 

(v) The expected oxygen (OviI and Ovill) line counts depend 
significantly on the choice of the feedback models. In particular, 
the number of detectable line counts is increased by a factor of 

3 by G Ws a nd b y 20 pe r cent when assuming a top-heavy IMF 
dArimoto & Yoshiiill987h . We obtain also a significant enhance- 
ment by increasing the mass resolution of the simulation. 

(vi) With our simulations we predict a number of detectable oxy- 
gen emission lines lower by a fac tor of 3-4 with respect to the pre- 
dictions of other works dCen & Fang 2006; Takei et al. 2011) due 
to differences in the implementation of chemical enrichment in the 
hydrodynamical codes. This has to be subject of further investiga- 
tions and has to be considered a part of current uncertainties in the 
modelisation of the WHIM properties. 

As a conclusion, our work confirms the importance of the 
use of hydrodynamical simulations in the modelisation of the LSS 
properties, together with the comparison with observations. In par- 
ticular, we highlighted how the new generation of X-ray telescopes 
could provide constraints on both the missing baryon component 
and the mechanisms of star-formation and feedback that affect the 
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WHIM. In this framework, it will be interesting to extend the re- 
sults of our work also to the absorption properties in the X-ray and 
UV frequency range. 
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